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In  probabilistic  design  of  materials  for  fracture-critical  components  in  modern  military  turbine  engines,  a 
typical  maximum  design  target  risk  (DTR)  is  5  x  10-8  component  failures/engine  flight  hour.  This  metric 
underscores  the  essential  role  of  safety  in  a  design  process  that  simultaneously  strives  to  achieve  perfor¬ 
mance,  efficiency,  reliability,  and  affordability  throughout  the  life  cycle  of  the  engine.  Traditionally,  the 
design  and  life  management  approaches  for  engine  materials  have  typically  relied  on  extensive  testing 
programs  to  produce  large  databases  of  fatigue  data,  from  which  statistically  based  life  limits  are  derived 
by  extrapolation  from  the  mean  fatigue  behavior.  However,  we  have  found  that  the  statistical  behavior  of 
fatigue  lifetimes  under  a  given  test  condition  often  exhibits  a  bimodal  form,  and  that  the  trends  in  mean 
vs.  minimum  fatigue  lifetime  typically  respond  differently  to  loading  or  to  microstructural  variables. 
Under  such  circumstances,  the  underlying  life-limiting  mechanisms  appear  to  exhibit  a  probabilistic 
microstructural  hierarchy  in  fatigue  resistance  that  is  controlled  by  susceptibility  of  local  microstructural 
neighborhoods  to  early  damage  and  the  growth  of  small  cracks.  These  findings  suggest  that  significant 
opportunities  may  exist  for  reductions  in  uncertainty  in  materials  life-cycle  prediction  and  management, 
if  such  hierarchies  can  be  understood  and  controlled.  This  paper  explores  the  potential  implications  of 
these  findings,  and  a  number  of  possible  approaches  are  suggested  for  incorporating  the  insights  of 
life-limiting  fatigue  into  methods  of  integrated  computational  materials  engineering  (ICME)  to  support 
optimized  life-cycle  design  of  materials  and  components  in  turbine  engines.  Benefits  of  this  approach 
appear  to  include  substantial  improvements  in  model  accuracy,  coupled  with  reduced  requirements 
for  materials  testing,  potentially  leading  to  a  significant  reduction  in  the  time  and  cost  to  develop,  vali¬ 
date,  transition,  and  implement  new,  more  fatigue  resistant  alloys. 
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1.  Introduction 

l  A.  A  bright  future  for  materials  design 

Looking  to  the  future  of  turbine  engine  materials,  it  appears  that 
there  are  many  potential  breakthroughs  in  the  underlying  science 
and  technology  of  materials  design  that  could  offer  important 
opportunities  for  improvement  beyond  today’s  engineering  state- 
of-the-art.  For  example,  future  designs  are  likely  to  benefit  from  a 
wave  of  anticipated  advancements  in  probabilistic  modeling  and 
simulation,  micro-experimentation,  high  fidelity  materials  charac¬ 
terization,  mechanism-based  modeling,  and  information  science 
and  technology,  all  of  which  are  areas  of  enhanced  investment 
under  a  number  of  major  new  initiatives  in  science  and  technology 
[1-4].  Through  unprecedented  levels  of  model  integration  among  a 
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range  of  science  and  engineering  disciplines,  it  should  be  possible  to 
design  future  propulsion  systems  to  be  more  adaptable  to  new 
usage  requirements,  to  have  higher  fuel  efficiency,  to  be  inherently 
more  reliable  and  verifiable,  to  be  manufactured  and  certified  more 
quickly,  and  to  require  far  less  logistics  support  over  the  system’s 
full  life  cycle. 

In  the  relatively  near  term,  continuous  recording  of  engine  flight 
and  performance  data  streams  is  expected  to  give  operators  and 
logisticians  unsurpassed  knowledge  of  the  usage  of  individual  en¬ 
gines.  Using  comprehensive  engine  performance  and  stress  analy¬ 
sis  models  of  the  engine,  one  can  then  calculate  location-specific 
stresses  and  temperatures  in  individual  engine  components  for 
use  in  life  analysis.  In  the  longer  term,  it  should  be  possible  to 
use  diagnostics  of  structural  and  materials  health  as  part  of  a  com¬ 
prehensive  philosophy  of  condition  based  management  +  progno¬ 
sis  (CBM+)  [5-9].  In  the  longer  term,  a  key  goal  is  to  pursue  the 
concept  of  aircraft  digital  twin  to  represent  the  evolving  health 
state  of  each  individual  asset  [10]. 
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Undergirding  these  forecasted  advancements,  the  burgeoning 
field  of  integrated  computational  materials  engineering  (ICME) 
[11,12]  is  expected  to  enable  game-changing  materials-design  con¬ 
cepts,  including  the  ability  to  tailor  location-specific  material  prop¬ 
erties  that  are  optimized  to  meet  component-specific  life-cycle 
design  requirements.  Perhaps  the  greatest  underlying  challenge 
for  the  science  and  engineering  community,  however,  will  be  to  ex¬ 
ploit  the  full  potential  of  a  material’s  capabilities,  while  simulta¬ 
neously  assuring  safety,  reliability,  and  affordability  over  the 
entire  system  life  cycle.  In  this  regard,  the  growing  importance  of 
affordability  was  recently  underscored  in  the  US  Air  Force  forecast 
report:  Technology  Horizons ,  Essential  Focus  Areas  for  Air  Force  S  and 
T  Investment,  which  noted  that  .  .Air  Force  science  and  technology 
over  the  next  decade  will  need  to  focus  as  much  on  advancing 
technologies  that  enable  reduced  Air  Force  operating  costs  as  on 
technologies  that  support  more  traditional  development  of  new 
systems  or  capabilities.”  [13]. 

To  help  meet  such  challenges,  the  engine  structural  design  com¬ 
munity  is  immersed  in  building,  improving,  and  exploiting  compre¬ 
hensive  probabilistic  design  tools  for  engine  components  (e.g., 
[14]).  Underpinning  such  computer  codes  are  the  variabilities  and 
uncertainties  in  the  input  parameters  and  data,  including  probabi¬ 
listic  descriptions  of  engine  missions/loads,  aerodynamic  flow, 
thermal  and  stress  analysis,  combustion  effects,  and  materials  per¬ 
formance  in  harsh  environments.  At  the  structures  and  system  le¬ 
vel,  new  concepts  for  robust  design  [15]  and  prognosis  [5]  are 
emerging.  In  addition,  recent  advancements  in  Bayesian  updating 
offers  a  rigorous  approach  to  refine  model  accuracy  based  on  con¬ 
tinuous  feedback  from  diagnostic,  maintenance,  and  operational 
data  [16].  While  it  may  take  decades  to  fully  develop  and  imple¬ 
ment  the  science  and  technology  needed  to  achieve  the  vision  out¬ 
lined  above,  there  are  many  incremental  steps  that  can  facilitate 
transition  and  which  offer  significant  nearer-term  benefits.  Perhaps 
the  key  initial  step  is  the  adoption  of  a  modular,  risk-based  model¬ 
ing  framework  in  which  input  data  and  assumptions  are  repre¬ 
sented  probabilistically.  Under  this  scenario,  each  input  to  the 
system-level  model  would  be  represented  by  one  or  more  sub¬ 
models,  which  could  be  individually  developed  and  validated.  Thus, 
the  set  of  models  of  the  usage,  structural  response,  and  materials 
performance  could  be  continually  improved,  as  the  underlying  sci¬ 
ence  advances  from  its  current  empirical  representation  to  a  more 


fundamental  physics  base.  The  intent  of  the  current  paper  is  to 
highlight  an  area  of  progress  in  understanding  uncertainty  in  the  fa¬ 
tigue  lifetime  of  a  selected  titanium  alloy  and  to  outline  related 
opportunities  in  the  area  of  integrated  computational  materials 
engineering. 

1.2.  Case  study:  turbine  engine  materials  -  building  on  a  history  of 
design  and  operation 

Life-cycle  management  procedures  used  for  today’s  advanced 
turbine  engines  have  achieved  a  superb  record  of  safety,  and  the 
underlying  models  and  methods  have  evolved  over  decades  of  de¬ 
sign  and  operation.  This  process  has  culminated  in  validated  engi¬ 
neering  standards  of  practice,  supported  by  large  databases  of 
material  properties.  Consequently,  these  engineering  practices 
are  highly  tolerant  to  variations  and  uncertainties  in  materials 
and  structures,  and  catastrophic  fatigue  failures  of  components 
have  been  virtually  eliminated.  At  the  same  time,  owing  largely 
to  the  mature  age  of  some  of  the  fleets,  engine  life-cycle  costs 
are  growing  at  an  annual  rate  that  is  well  above  inflation,  and  this 
cost  burden  threatens  to  stifle  efforts  to  recapitalize  and  modernize 
future  engine  fleets  [17].  Given  such  affordability  challenges,  there 
is  a  growing  need  to  reduce  life-cycle  costs  for  both  legacy  and  fu¬ 
ture  engines.  However,  history  indicates  that  as  much  as  80%  of  an 
engine’s  life  cycle  cost  is  established  during  the  design  and  early 
development  of  the  engine,  and  the  cost  to  correct  design  defects 
afterward  becomes  exponentially  greater  as  development  and 
deployment  proceed  [17]. 

Within  the  US  Department  of  Defense,  the  decades  of  engine  de¬ 
sign  and  engineering  experience  have  molded  the  specifications 
known  as  the  Propulsion  System  Integrity  Program  (PSIP)  [18] 
and  its  associated  document  [19],  which  guide  the  life-cycle  design 
and  management  of  turbine  engines.  The  process  of  PSIP  design, 
depicted  schematically  in  Fig.  1,  is  the  beginning  of  a  comprehen¬ 
sive  life-cycle  design  and  logistics  management  process  that  in¬ 
cludes:  (i)  design  information  and  development  planning;  (ii) 
design  characterization  and  analysis;  (iii)  components  and  subsys¬ 
tems  testing;  (iv)  propulsion  system  ground  and  flight  testing,  and 
(v)  the  engine  life  management  plan.  Engine-material  develop¬ 
ment,  selection,  characterization,  modeling,  and  validation  are  con¬ 
centrated  in  tasks  (i)  and  (ii),  and  first-tier  material  properties, 
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such  as  strength,  stiffness,  and  density,  tend  to  drive  initial  design. 
Second-tier  material  properties,  such  as  cycle-,  time-,  temperature- 
,  and  environmental-dependence,  are  prominent  factors  through¬ 
out  the  design  and  life  cycle,  and  these  properties  typically  define 
engine  overhaul  intervals,  control  the  average  time  on  wing,  and 
dictate  the  ultimate  lifetimes  of  components. 

For  safety-  or  mission-critical  components,  such  as  rotors  and 
rotor  spacers,  the  PSIP  design  life  is  driven  primarily  by  major  load 
cycles  produced  by  variations  in  engine  rotational  speed,  wherein  a 
ground-air-ground  flight  is  counted  as  one  cycle,  and  minor  throt¬ 
tle  excursions  during  flight  count  as  partial  cycles.  Such  cycles  may 
be  summarized  as  total  accumulated  cycles  (TACs),  and  a  typical 
flight  produces  approximately  2-3  TACs/h.  To  satisfy  performance 
and  safety  requirements,  a  design  life  for  major  rotating  compo¬ 
nents  in  high  performance  engines  is  nominally  8000  TACs,  with 
a  complete  overhaul  of  the  engine  prescribed  at  4000  TACs.  Given 
these  target  lifetimes,  the  associated  maximum  stresses  in  rotor 
components  tend  to  approach  the  material  yield  strength  at  key 
locations,  such  as  boltholes  and  other  stress  concentrators,  and  rel¬ 
atively  high  maximum  stresses  occur  throughout  the  components. 

With  regard  to  the  materials  life  cycle,  it  is  well  known  that  ad¬ 
vanced  Ti  and  Ni  alloys  in  turbine  engine  rotors  exhibit  a  substan¬ 
tial  variation  in  their  long-term  cyclic-  and  time-dependent 
properties,  such  as  resistance  to  fatigue  damage  and,  to  a  lesser  ex¬ 
tent,  resistance  to  creep  deformation.  The  PSIP  specifications  deal 
with  such  variation  by  managing  to  a  criterion  based  on  statistical 
lower-bound  behavior.  For  example,  for  design  of  safety-  or  mis¬ 
sion-critical  components,  these  specifications  impose  simulta¬ 
neous  requirements  for  (i)  a  statistical  minimum  safe-life  fatigue 
and  (ii)  damage  tolerance  to  small  surface  or  corner  cracks.  Under 
the  safe-life  criterion,  the  minimum  LCF  life  is  defined  as  -3  stan¬ 
dard  deviations  from  the  mean  value  of  the  material’s  low  cycle  fa¬ 
tigue  lifetime,  or  alternatively,  the  1/1000  (B0.1)  lower-bound 
capability,  based  on  a  database  of  tests  over  the  range  of  engine  de¬ 
sign  conditions.  Thus,  the  PSIP  specification  virtually  guarantees 
material  safety  under  low  cycle  fatigue,  although  the  reliance  on 
large  statistical  databases  can  be  costly  and  time  consuming. 

Engines  are  designed  to  meet  targeted  mission  profiles  of  loads, 
but  actual  mission  usage  invariably  changes  as  the  operational 
needs  change,  making  tracking  of  actual  mission  usage  a  key  input 
to  improved  life  analysis.  Future  engines  are  to  be  designed  with 
sophisticated  usage-tracking  capabilities,  and  the  stresses  and 
temperatures  of  engine  components  are,  therefore,  expected  to 
be  much  more  well  defined  than  is  possible  today.  However,  the 
understanding  and  control  of  variability  in  the  material’s  response 
to  such  loading  is  much  less  mature,  and  represents  a  significant, 
although  challenging,  opportunity  for  improvement. 

In  the  quest  for  improved  models  to  address  the  goals  of  safety, 
reliability,  and  affordability,  safety  will  certainly  remain  as  the  first 
requirement  for  manned  aircraft.  For  example,  as  noted  in  the  PSIP 
specification,  the  limit  for  the  design  target  risk  (DTR)  for  probabi¬ 
listic  design  of  a  component  “was  set  at  a  value  of  5  x  10-8  pre¬ 
dicted  occurrences  per  EFH  [Engine  Flight  Hour]  for  a  recent, 
advanced  single-engine  fighter.”  Similarly,  specifications  for  com¬ 
mercial  engine  components  set  the  maximum  DTR  at  1  x  1 0-9/ 
flight  [20].  Such  rigorous  safety  requirements  underscore  the 
importance  of  fatigue  life  limits  in  design  and  operation  of  turbine 
engines,  and  highlight  the  challenges  facing  integrated  computa¬ 
tional  materials  engineering  for  engine  alloys. 

1.3.  Life  limits  and  competing-mode  of  fatigue 

Under  previous  efforts  on  life-limiting  material  behavior,  we 
have  found  that  the  variability  in  fatigue  response  of  engine  alloys 
at  a  given  stress  level  often  tends  to  separate  into  bimodal  distribu¬ 
tions  representing  minimum  behavior  and  a  mean-life  dominated 


behavior,  termed  Type  I  and  Type  II  behavior,  respectively.  Thus, 
the  probability  density  representing  a  material’s  safe-life  fatigue 
capability,  depicted  in  Fig.  1,  was  often  composed  of  two 
underlying  probability  densities,  as  shown  schematically  in 
Fig.  2.  This  separation  in  fatigue  lifetimes,  known  as  bimodal  or 
competing-modes  fatigue,  has  been  observed  in  a  wide  range  of  al¬ 
loys,  including  the  superalloys:  Rene’95  [21,22],  Rene’88DT  [23], 
INI 00  [8,24-29],  Waspaloy  [30],  the  single  crystal  alloy  PWA 
1484  [31],  the  titanium  alloys  Ti-10-2-3  [23,32-34],  Ti-6Al-2Sn- 
4Zr-6Mo  [35-42],  Ti-6A1-4V  [43-47],  gamma  titanium  aluminides 
[48,49],  the  aluminum  alloy  7075-T651  [50],  and  others.  Although 
such  a  separation  of  fatigue  response  has  been  known  for  some 
time  [51  ],  the  significance  of  this  behavior  has  not  yet  been  gener¬ 
ally  captured  in  the  strategies  for  fatigue  design  of  turbine  engine 
materials. 

In  exploring  this  behavior  in  turbine  engine  alloys,  we  have 
developed  a  physically-based  approach  for  describing  fatigue  var¬ 
iability,  and  this  approach  has  been  integrated  into  a  Monte  Carlo 
probabilistic  life  prediction  model  for  materials  [35,40,52].  The  pri¬ 
mary  tenets  of  our  fatigue  variability  description  are  summarized 
by  [40]:  (i)  under  nominal  microstructural  and  loading  conditions, 
a  hierarchy  of  local  deformation  heterogeneities  develop  in  a  fati¬ 
gue  sample,  corresponding  to  grain-scale  cyclic  micro-plasticity  in 
certain  microstructural  arrangements  or  features;  (ii)  as  a  result  of 
this  hierarchy,  a  probability  exists  of  an  extreme  microstructural 
arrangement  that  may  initiate  a  crack-growth  dominated  failure, 
producing  a  life-limiting  failure  distribution  for  the  material;  and 
(iii)  there  are  often  different  degrees  of  influence  of  a  given  variable 
on  minimum  vs.  the  mean  lifetimes.  For  example,  the  effect  of 
stress  level  on  mean  vs.  life-limiting  fatigue  typically  produces 
quite  different  responses,  and  such  differing  response  of  minimum 
and  mean  fatigue  behavior  has  been  observed  for  other  important 
variables,  including  microstructure,  surface  treatment,  tempera¬ 
ture,  and  dwell  time.  Furthermore,  this  behavior  appears  to  mani¬ 
fest  itself  in  the  presence  of  notches  [53].  This  understanding, 
when  applied  to  fatigue  design  of  materials,  suggests  that 
compositional,  processing,  surface  treatment,  and  microstructural 
modifications  should  consider  the  lower-lifetime  limit,  and  that 
the  life-limiting  behavior  may  not  be  accurately  represented  by 
purely  statistical,  mean-fatigue-based,  approaches. 

The  goal  of  the  current  paper  is  to  (i)  overview  general  findings 
on  fatigue  life  limits  of  alloys  used  in  turbine-engine  rotor 
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Fig.  2.  Expanding  on  the  concept  of  “Expected  Material  Properties”  highlighted  at 
the  upper  left  of  Fig.  1,  it  is  commonly  assumed  that  fatigue  variability  is  well 
described  by  a  single  probability  distribution  (light  blue),  as  show  above. 
Traditional  methods  of  design  and  life  management  derive  material  life  limits  by 
extrapolation  to  the  lower  tail  of  the  single  distribution.  However,  it  now  appears 
that  such  fatigue  lifetime  behavior  may  actually  arise  from  at  least  two  separate 
underlying  probability  distributions.  In  this  sense,  the  identification  of  two 
distributions  appears  to  offer  substantial  opportunity  for  improved  understanding 
of  the  fundamental  material  behavior  and  a  reduction  in  predictive  uncertainty. 
(For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 
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components,  (ii)  to  highlight  the  effects  and  implications  of  such 
fatigue  life  limits  for  a  specific  high  strength  titanium  alloy,  and 
(iii)  to  identify  possible  avenues  for  improvement  in  materials- 
model  accuracy.  Thus,  the  focus  will  be  on  challenges  and  opportu¬ 
nities  to  reduce  variability  in  performance,  and  uncertainty  of  pre¬ 
dictions,  of  fatigue  life-cycle  capability  of  turbine-engine  rotor 
materials,  consistent  with  the  demands  of  a  high-performance  sin¬ 
gle-engine  aircraft. 

2.  Materials  and  experimental  methods 

The  focus  of  this  investigation  was  a  pancake  forging  of  alloy  Ti- 
6Al-2Sn-4Zr-6Mo(Ti-6246)  processed  and  heat-treated  to  a  du¬ 
plex  a  +  p  phase  microstructural  condition,  as  shown  in  Fig.  3. 
The  microstructure  consisted  of  primary,  equiaxed  a  particles 
(ap),  and  lamellar  colonies  of  fine  a//?  laths,  which  transformed 
from  prior- 13  grains  during  heat  treatment.  The  average  size  of 
the  ap  was  about  4  pm  and  that  of  the  colonies  was  about 
15  pm.  The  volume  fractions  of  ap  and  the  colony  phases  were 
approximately  0.3  and  0.7,  respectively.  This  produced  a  bulk  yield 
strength  of  1140  MPa  and  an  ultimate  strength  of  1250  MPa. 

Closely  spaced,  cylindrical  specimen  blanks  of  this  material 
were  machined  from  the  forging,  and  sections  of  Ti-6A1-4V  rod 
were  inertia  welded  to  the  ends  as  extensions  for  gripping.  Subse¬ 
quently,  the  inertia-welded  blanks  were  machined  and  finished  by 
low-stress  grinding  into  either  cylindrical  dog-bone  specimens  for 
fatigue  testing,  or  similar  small-crack  specimens,  which  were  fab¬ 
ricated  with  surface  flats  to  accommodate  crack  growth  measure¬ 
ments  using  standard  acetate  replication  procedures.  Surface 
residual  stresses  from  machining  were  removed  by  electropolish¬ 
ing  the  machined  specimens  to  remove  a  50  pm  thick  surface  layer 
on  the  specimen  gage  sections.  The  small-crack  tests  were  per¬ 
formed  using  either  natural  initiation  or  by  initiation  from  multiple 
small  surface  notches,  which  had  been  machined  into  the  surface 
using  a  focused  ion  beam  (FIB)  to  create  residual-stress-free  crack 
starters. 

To  assess  the  inherent  effect  of  the  Ti-6246  alloy  microstructure 
on  the  variability  in  fatigue  response,  a  series  of  tension-tension 
fatigue  tests  was  performed  at  each  of  nine  stress  levels  under  a 
stress  ratio  of  minimum/maximum  stress  of  0.05.  The  specific 
maximum  stresses  for  these  tests  were  1040,  925,  900,  860,  820, 
700,  650,  600,  and  550  MPa.  The  test  frequency  was  20  Hz  for  max¬ 
imum  stresses  above  750  MPa,  while  tests  below  this  stress  were 


20pm 


Fig.  3.  Nominal  duplex  microstructure  of  Ti-6Al-2Sn-4Zr-6Mo  alloy  processed  to 
in  the  oc  +  f  phase  field,  containing  approximately  30%  primary  a  grains  and  fine  a- 
phase  platelets  in  a  matrix  of  transformed  f  phase. 


performed  in  ultrasonic  resonance  at  20  kHz.  At  the  highest  stres¬ 
ses,  surface-connected  a  particles  typically  served  as  the  primary 
sites  for  crack  initiation  (e.g.,  Fig.  4),  while  some  of  the  specimens 
failed  from  internal  initiation  sites  under  maximum  stresses  of 
860  MPa  or  lower.  Further  details  on  the  microstructure,  experi¬ 
mental  methods,  results,  and  analysis  are  provided  in  [37,39- 
41,54],  which  served  as  the  backbone  of  the  results  to  be  presented 
below. 

3.  Results 

Fig.  5  presents  an  S-N  plot  of  the  resulting  fatigue  lifetimes  for 
all  stress  levels  (crmax  of  550-1040  MPa)  and  both  frequencies.  As 
shown,  the  variability  in  fatigue  lifetime  exceeded  three  orders  of 
magnitude  at  the  intermediate  stress  levels,  while  the  variability 
in  lifetime  at  the  highest  stress  level  (1040  MPa)  was  extremely 
limited.  The  increase  in  variability  with  reducing  stress  has  been 
observed  widely  in  the  literature,  although  the  extent  of  the  in¬ 
crease  is  likely  maximized  for  the  current  tests,  due  to  the  use  of 


10pm 

Fig.  4.  Typical  fatigue  crack  origin  (lower  edge)  showing  a  single,  surface- 
connected,  alpha  grain  facet  where  crack  nucleation  occurred. 


Fig.  5.  Ti-6246  stress  vs.  life  (S-N)  with  plot  showing  individual  specimen  lifetimes. 
Tests  above  750  MPa  performed  at  20  Hz,  and  tests  below  750  MPa  performed  at 
20  kHz. 
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residual-stress-free  electropolished  surfaces.  This  surface  condi¬ 
tion  was  chosen,  however,  precisely  because  it  eliminated  effects 
of  surface  irregularities  and  residual  stresses  and  helped  to  reveal 
the  inherent  effects  of  alloy  microstructure.  The  nature  of  the  sta¬ 
tistical  variations  in  the  fatigue  lifetimes  at  each  stress  level  is  elu¬ 
cidated  in  the  probability  plot  of  Fig.  6.  As  shown,  the  lifetimes 
from  the  highest  stresses  exhibited  a  straight-line  trend  on  the  log¬ 
normal  probability  scale,  while  bimodal  trends  of  the  data  tended 
to  emerge  at  intermediate  stress  levels.  As  shown  on  an  expanded 
scale  in  Fig.  7,  the  bimodal  behavior  was  most  evident  for  the  fati¬ 
gue  datasets  at  crmax  =  860  and  820  MPa.  For  each  of  these  datasets, 
three  regions  of  fatigue  behavior  are  distinct  (Type  I,  transition, 
and  Type  II),  with  the  ratio  of  the  data  points  in  the  Type  I  vs.  Type 
II  regions  being  stress-dependent.  The  transition  points  appear  to 
represent  the  short-life  region  of  the  Type  II  behavior,  indicating 
a  substantial  crack  initiation  lifetime.  However,  no  standard  prob¬ 
ability  density(s)  adequately  described  the  transition  region,  sug¬ 
gesting  that  the  transition  may  represent  a  change  in  initiation 
mechanism.  More  work  is  required  to  clarify  this  speculation, 
however. 

These  separate  regions  of  fatigue  response  were  used  as  a  basis 
for  isolating  the  short-  and  long-life  trends  alone  (eliminating  the 
transition  regions),  as  re-plotted  in  Fig.  8.  As  shown,  when  the  Type 
I  and  Type  II  regions  are  plotted  separately,  the  data  appear  nom¬ 
inally  straight  on  the  lognormal  probability  scale,  indicating  that 


Fig.  6.  Lognormal  probability  plot  of  fatigue  lifetime  data  for  individual  stress  levels 
shown  in  Fig.  5.  From  left  to  right:  <rmax=  1040,  925,  900,  860,  820,  700,  650,  600, 
and  550  MPa. 


Fig.  7.  Lognormal  probability  plot  showing  820  MPa  and  860  MPa  data  from  Fig.  5, 
highlighting  three  regions  of  the  bimodal  trend  (Type  I,  transition,  Type  II). 


each  region  is  unimodal.  Note  that  the  data  from  the  crmax  =  860 
and  820  MPa  conditions  are  plotted  twice,  since  both  Type  I  and 
Type  II  regions  were  clearly  present  at  these  stress  levels.  In  com¬ 
parison  with  the  full  range  of  fatigue  variability  evident  in  the 
complete  datasets  of  Fig.  6,  Fig.  8  shows  that  the  two  regions  of 
behavior  tend  to  follow  distinctly  different  trends,  appearing  to 
represent  a  life-limiting  (Type  I)  behavior,  plus  a  final  period  of 
fatigue  exhaustion  where  all  specimens  eventually  failed.  Interest¬ 
ingly,  in  terms  of  the  lognormal  probability  scale  of  the  plot,  there 
was  limited  variation  of  the  slopes  of  the  data  at  the  individual 
stress  levels  in  the  two  regions  of  behavior  (life-limiting  and  fati¬ 
gue-exhaustion),  and  the  general  slopes  of  the  two  regions  were 
remarkably  similar.  Although  all  specimens  failed  (no  run-out 
specimens),  it  appears  that  the  longest-lived  specimens  may  have 
been  approaching  a  fatigue  run-out  of  infinite  lifetime. 

Fig.  9  shows  the  data  from  the  820  MPa  tests,  and  includes  a 
lognormal  probability  regression  line  with  95%  confidence  inter¬ 
vals  on  the  regression.  Clearly,  such  a  regression  is  inappropriate 
for  bimodal  data.  However,  materials  and  component  designers 
rarely,  if  ever,  have  the  luxury  of  having  so  many  repeated  tests 
at  an  individual  fatigue  condition.  With  fewer  test  results,  the  bi¬ 
modal  behavior  would  probably  not  be  evident,  and  a  regression 


Cycles  to  Failure 

Fig.  8.  Lognormal  probability  plot  of  Type  I  and  Type  II  fatigue  lifetime  data  alone 
for  individual  stress  levels  shown  in  Fig.  5.  Note  that  820  and  860  MPa  data  each 
appear  twice,  since  these  stress  levels  exhibited  both  Type  I  and  Type  II  behavior.  As 
shown  in  Fig.  6,  data  from  700  and  900  MPa  tests  also  exhibited  both  types  of 
behavior,  but  to  a  lesser  extent  than  at  820  and  860  MPa. 


Fig.  9.  Lognormal  probability  plot  of  lifetime  data  for  tests  with  <7max  =  820  MPa, 
showing  a  regression  extrapolation  to  the  B0.1  (1/1000)  probability  of  failure,  along 
with  95%  confidence  bands  for  the  fit. 
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similar  to  the  one  shown  would  likely  appear  reasonable.  In  this 
case,  one  can  see  that  extrapolation  to  the  B0.1  probability  of  fail¬ 
ure  would  give  a  lifetime  of  approximately  3500  cycles,  while  the 
corresponding  lower  bound  of  the  95%  confidence  interval  was 
about  721  cycles.  Alternatively,  taking  maximum  benefit  of  the 
large  number  of  data  points  at  this  condition  allows  one  to  focus 
on  the  single,  short-life  (Type  I)  distribution.  Fig.  10  shows  a  prob¬ 
ability  plot  of  the  6  shortest-life  data  alone,  and  in  this  instance  the 
estimated  B0.1  life  was  greater  than  104  cycles,  and  the  minimum 
95%  confidence  estimate  was  approximately  5660  cycles,  or  about 
8  times  the  corresponding  value  indicated  by  the  unimodal  analy¬ 
sis  of  the  full  data  set.  Thus,  this  example  implies  that  existence  of 
bimodal  fatigue  behavior  may  indicate  that  the  material’s  inherent 
minimum  lifetime  is  substantially  greater  than  otherwise  thought. 


4.  Analysis 

4A.  Modeling  bimodal  fatigue 

A  key  opportunity  in  the  design  of  future  materials  is  to  under¬ 
stand,  isolate,  and  model  the  sources  of  variability  in  fatigue  dam¬ 
age  and  lifetime.  As  noted  earlier,  we  previously  established  a 
general  computational  framework  for  physics-based  understand¬ 
ing  of  fatigue  variability  [24,40].  In  this  approach,  the  lifetime  dis¬ 
tribution  was  modeled  as  a  superposition  of  the  probability 
densities  of  the  crack-growth-dominated  lifetime  (Type  I)  and 
the  crack-initiation-dominated,  mean  lifetime  (Type  II).  The  under¬ 
lying  initiation  and  cracking  mechanisms  for  this  material  have 
been  reported  elsewhere,  and  a  key  factor  in  the  fatigue  variability 
of  this  material  was  microtexture  [40-42,54-56].  The  total  lifetime 
density  was  represented  by: 

ft(N)=P]fm+PJm(N)  (1) 

where  /t(N),  f  (N),  and  /m(N)  are  probability  densities  representing 
the  total  lifetime,  crack  growth  lifetime,  and  the  mean,  initiation- 
dominated  lifetime,  respectively.  The  factors,  pi  and  pm  are  the  prob¬ 
abilities  of  occurrence  of  the  respective  responses.  Salient  features  of 
the  model  are  (i)  that  the  effect  of  any  variable  on  lifetime  distribu¬ 
tion  is  separated  into  its  influence  on  small-plus-large  crack  growth 
vs.  the  mean  lifetime  behaviors,  and  (ii)  that  the  lifetime  distribution 
is  modeled  as  a  superposition  of  these  two  responses.  From  a  design 
life  perspective  (for  example,  a  1/1000,  B0.1,  lifetime),  it  can  be 
shown  that  the  prediction  of  life  limit  by  the  crack  growth  lifetime 


density, /i(N)  (i.e.,  in  the  limit  of  pi  =  1 )  forms  a  lower  bound  of  pre¬ 
dictions  by  the  bimodal  model  [35],  as  discussed  below. 

It  is  well  known  that  the  total  fatigue  life  of  a  specimen  is  com¬ 
posed  of  the  sum  of  the  cycles  from  the  successive  stages  of  dam¬ 
age  initiation  and  crack  growth  to  failure,  such  that  the  total 
fatigue  life  (NT)  is  given  by: 

Nt^Nj+Nsc  +  Nlc  (2) 


where  Nx,  Nsc,  and  NLC  are  the  cycles  spent  in  crack  initiation,  small 
crack  growth,  and  large-crack  growth,  respectively.  It  is  also  recog¬ 
nized  that  N[  may  be  further  subdivided  into  a  crack-incubation  plus 
a  crack-nucleation  phase  in  some  instances,  and  that  small  cracks 
may  be  defined  as  being  microstructurally,  physically,  mechani¬ 
cally,  and/or  chemically  small  [57-59].  For  a  given  material  and 
loading  condition,  the  variability  in  total  fatigue  lifetime,  VT,  is  a 
sum  of  the  variabilities  of  the  individual  stages  of  damage  and 
cracking,  and  variability  of  the  individual  terms  in  Eq.  (2)  tends  to 
reduce  from  left  to  right,  such  that  V{NX)  >  V{NSc)  >  V{NLC). 

As  noted  earlier,  for  a  typical  8000-cycle  design  life  of  turbine- 
engine  rotor  materials,  the  corresponding  maximum  stresses  ap¬ 
proach  the  material’s  yield  stress.  Under  these  circumstances,  it 
appears  that  the  number  of  cycles  to  crack  initiation  from  a  micro- 
structural-scale  feature  approaches  zero  for  the  life-limiting  fati¬ 
gue.  Assuming  that  Nx  ~  0,  one  can  calculate  the  expected 
lifetime  by  integrating  the  equation: 
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where  the  two  integrals  represent  the  life  spent  in  small-crack  and 
large-crack  growth,  respectively,  ax  is  the  crack  initiation  depth 
from  a  microstructural  or  mechanical  defect,  aSc  is  the  transition 
crack  size  between  small-  and  large-crack  growth,  and  /sc  and  /LC 
are  the  crack  growth  functions  for  the  small-  and  large-crack  size 
ranges,  respectively.  For  the  Monte  Carlo  model  of  the  life-limiting 
fatigue  process,  the  distribution  of  crack  initiation  sizes  {ax)  was 
represented  by  the  red  curve  of  Fig.  11.  For  the  data  presented,  ax 
was  determined  from  a  combination  of  surface  measurements  via 
standard  acetate  replica  methods,  combined  with  fractography  of 
the  crack  initiating  features,  which  were  typically  individual  equi- 
axed  primary  a-phase  particles  ~4  um  in  diameter,  such  as  the 
one  shown  in  Fig.  4.  Details  on  the  initiation  fractography  are  avail¬ 
able  in  [39-41,52,55].  Analysis  of  the  crack  initiation  sites  revealed 
that  the  small  cracks  appeared  to  assume  an  approximately 
semicircular  shape  very  quickly.  As  such,  ax  was  defined  as  the 


Fig.  10.  Lognormal  regression  of  the  Type  I,  shortest-life,  points  alone,  showing  that 
the  slope  of  the  regression  is  much  steeper  and  the  B0.1  extrapolated  life  is  much 
greater  than  shown  from  the  full  dataset  in  Fig.  8. 


Fig.  11.  Histograms  of  primary  a  particle  sizes  (ap)  sizes  (blue)  as  measured  from 
metallurgical  section  and  sizes  of  crack  initiation  regions  (red)  measured  as 
projected  normal  to  the  fracture  surface.  (For  interpretation  of  the  references  to 
colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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equivalent  radius  of  the  initiating  feature,  and  the  crack  was  as¬ 
sumed  to  be  semicircular,  since  correlation  with  specific  shapes  of 
initiation  particles  was  beyond  the  scope  of  the  current  project. 
General  information  on  shape  effects  for  small  cracks  in  other  mate¬ 
rials  is  available  in  [60-64].  A  statistical  representation  of  the  small- 
crack  growth  behavior  shown  in  Fig.  12  (ft  =  0.05)  was  obtained 
using  maximum  likelihood  methods  to  obtain  a  correlated  Paris- 
law  function  to  represent  the  variability  in  crack  growth  behavior, 
while  the  large-crack  behavior  was  treated  as  deterministic  (mean 
value).  Since  the  small  crack  data  were  acquired  at  the  same  stress 
ratio  as  the  smooth-bar  fatigue  specimens,  no  correction  for  ft  was 
used.  Moreover,  since  the  small  crack  data  exhibited  no  apparent 
threshold  for  crack  growth  (A/Cth),  this  value  was  assumed  to  be 
zero  in  the  model  calculations.  For  each  stress  level,  at  least 
10,000  Monte  Carlo  simulations  were  used  to  calculate  a  distribu¬ 
tion  of  expected  minimum  total  lifetimes,  assuming  N[  =  0.  The 
Monte  Carlo  simulation  was  implemented  using  Mathematica®, 
and  the  stress  intensity  factor  calculation  was  based  on  the  solution 
for  a  surface  crack  in  a  round  bar  given  by  Forman  and  Shivakumar 
[65].  The  calculations  assumed  a  single  dominant  crack,  consistent 
with  the  experimental  findings,  and  the  failure  criterion  was  based 
on  the  effective  fracture  toughness.  Fig.  13  plots  these  predicted 
minimum  lifetimes  in  comparison  with  the  <rmax  =  820  MPa  life- 
limit  data  of  Fig.  10.  As  shown,  the  probabilistic  prediction  agrees 
extremely  well  with  the  life-limit  data,  which  tends  to  corroborate 
the  assumption  that  Nj  was  near  zero  for  these  tests. 

Fig.  14  re-plots  the  S-N  data  of  Fig.  5,  and  shows  a  schematic  of 
bimodal  regressions  of  the  fatigue  behavior  for  each  stress  level. 
The  three  lines  at  the  left  of  the  plot  represent  Monte  Carlo  model 
predictions  of  the  0.1%,  10%,  and  50%  probabilities  of  failure, 
respectively,  from  left  to  right.  As  shown,  all  of  the  fatigue  data 
are  bounded  approximately  by  the  P( failure)  =  10%  line.  In  addi¬ 
tion,  the  slopes  of  the  Monte  Carlo  predictions  vs.  stress  level  agree 
well  with  the  slope  of  the  Type  I  fatigue  lifetime  line  representing 
the  behavior  of  the  life-limiting  populations  alone,  as  defined  by 
the  linear  behavior  on  a  probability  plot.  Among  the  implications 
of  the  steeper  slope  for  life-limiting  fatigue  is  a  reduced  sensitivity 
to  stress  level,  as  compared  to  a  traditional  view  of  the  mean  S-N 
fatigue  behavior.  Additional  work  is  planned  to  explore,  validate, 
and  assess  the  generality  of  this  apparent  implication.  This  re¬ 
search  includes  work  on  micro-scale  fatigue  tests  to  help  charac¬ 
terize  modes  of  crystal  plasticity  leading  to  fatigue  damage  in  the 
Ti-6246  material  [66]. 

Also  as  shown  in  the  figure,  the  proportions  of  data  falling  in  the 
short-  (Type  I)  vs.  long-life  (Type  II)  regions  of  S-N  behavior  shift 


Fig.  12.  Plot  of  da/dN  vs.  A I<  for  several  naturally  initiated  small  surface  cracks 
measured  via  acetate  replication  of  the  specimen  surface,  and  regression  trend  from 
similar  data  on  large  cracks  in  conventional  compact  type  C(T)  specimens. 


Fig.  13.  Probability  plot  of  Monte  Carlo  model  simulation  of  Type  I  fatigue  lifetimes 
based  only  on  microstructure  (Fig.  12)  and  small  +  large  crack  (Fig.  13)  results 
plotted  vs.  the  short  life  data  from  the  <7max  =  820  MPa  experiment  (Fig.  9). 


Fig.  14.  S-N  plot  showing  (i)  individual  regression  lines  (blue  and  dashed)  to  Type  I 
and  Type  II  data,  (ii)  schematic  bimodal  probability  distribution  transitioning  from 
Type  I  to  Type  II  dominance  with  decreasing  stress,  and  (iii)  Monte  Carlo  predictions 
of  50%,  10%,  and  0.1%  probabilities  of  failure  for  the  minimum  lifetimes.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 

toward  longer  lifetimes  as  stress  is  reduced.  Recalling  Eq.  (2),  this 
is  equivalent  to  reducing  the  probability  (pi)  of  a  given  fatigue  data 
point  being  in  the  short  life  region.  Based  on  a  nonlinear  bimodal 
regression  of  the  fatigue  data  at  each  stress  level,  an  estimate  of 
the  dependence  ofpi  on  stress  level  is  plotted  in  Fig.  15.  As  shown, 
the  maximum  stress  at  which  the  short-  and  long-life  behaviors 
become  equally  likely  is  approximately  820  MPa,  which  is  about 
79%  of  the  yield  strength  for  the  Ti-6246  alloy  tested. 

Fig.  16a  presents  the  results  of  1000  Monte  Carlo  simulations  of 
crack  growth  behavior  at  820  MPa,  and  Fig.  16b  shows  an  ex¬ 
panded  region  of  the  same  plot.  As  shown,  for  the  distribution  of 
crack  initiation  depths  (up  to  approximately  15  pm  in  this  case), 
the  minimum  lifetime  is  more  than  14,000  cycles,  which  is  well 
in  excess  of  the  typical  8000  cycle  (TAC)  design  lifetime.  Recogniz¬ 
ing  that  the  estimate  of  p\  at  this  stress  level  is  0.5,  half  of  all  spec¬ 
imens  tested  at  this  condition  would  be  expected  to  lie  within  the 
short-life  distribution.  Since  traditional  low-cycle-fatigue  safe-life 
design  targets  the  1/1000  statistical  minimum  for  initiation  of  a 
0.40  mm  deep  semicircular  crack,  this  should  correspond  to  the 
minimum  values  of  the  Monte  Carlo  predictions.  As  such,  an  alter¬ 
native  interpretation  of  the  1/1000,  or  -3  standard  deviation,  LCF 
safe-life  criterion  is  that  for  these  cases,  in  all  likelihood,  beginning 
at  the  size  of  a  microstructural  feature  (~5-15  pm),  a  crack  would 
have  been  growing  for  most  of  the  safe-life  period.  There  is  a  nota¬ 
ble  caveat  in  this  regard,  however,  since  the  current  results  are  for 
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Fig.  15.  Probability  of  occurrence  of  Type  I  (short-life)  fatigue  based  on  bimodal 
regression  analysis  of  fatigue  lifetime  data  from  Fig  14. 


Fig.  16.  (a)  1000  Monte  Carlo  model  predictions  of  crack  depth  vs.  life  (a-N)  for  a 
maximum  stress  of  820  MPa,  and  (b)  and  the  same  predictions  on  an  expanded 
scale  that  shows  the  area  of  primary  interest  for  design  and  life  management  of 
rotor  components  designed  to  a  total  lifetime  of  8000  cycles  (TACs). 

a  stress-free  surface  condition,  while  actual  components  typically 
have  highly  compressive  surface  residual  stresses  from  machining 
plus  shot  peening,  as  discussed  below. 

4.2.  Effects  of  surface  residual  stresses  produced  by  shot  peening 

Surface  treatments  such  as  shot  peening  are  used  widely  to  re¬ 
tard  crack  initiation  and  propagation  in  the  surface  region,  and 
thereby  increase  the  fatigue  life.  However,  our  findings  from  other 
research  suggest  that  the  effect  of  shot  peening  on  the  lifetime  dis¬ 
tribution  should  be  represented  in  terms  of  its  effect  on  the  mean 
and  the  life-limiting  behavior  [42].  In  work  on  the  same  Ti-6246  as 
reported  in  the  current  paper,  shot-peening  increased  both  the 
minimum  and  mean  lifetimes  by  suppressing  crack-initiation  in 


the  surface  region  and  forcing  crack-initiation  to  the  subsurface 
in  many  samples  [42].  Although  failure  by  the  subsurface  mecha¬ 
nism  improved  the  mean  response,  the  life-limiting  failures  under 
shot-peening  occurred  by  a  different  surface-crack-initiation 
mechanism  in  which  cracking  was  found  to  originate  as  an  interac¬ 
tion  between  shot-peening  and  the  local  microstructure,  which 
created  local  regions  of  surface  micro-cracking.  These  surface  crack 
initiation  sizes  appeared  to  be  significantly  larger  than  the  average 
primary  a  grain  size  (~4  pm)  that  controlled  crack  initiation  in  the 
residual  stress  free  condition.  The  limiting  lifetimes  under  shot 
peening  were  only  about  1.7  times  longer  than  those  under  the 
residual-stress  free  condition.  In  other  words,  a  portion  of  the  ben¬ 
efit  from  compressive  surface  residual  stresses  induced  via  shot 
peening  was  partially  offset  by  the  creation  of  surface  cracks  from 
shot  peening  damage.  The  primary  benefit  of  the  residual  stress  ap¬ 
peared  to  be  due  to  the  stresses  acting  perpendicular  to  the  crack 
plane,  which  reduced  the  effective  range  of  crack  opening  and 
stress  intensity  factor  (AJQff). 

4.3.  Opportunities  for  improvement:  sources  of  variability  and 
uncertainty  in  material  life  limits 

Given  that  life  limits  play  such  a  crucial  role  in  the  design  and 
management  of  turbine  engine  rotor  materials,  it  appears  that 
the  results  outlined  in  this  paper  represent  a  mere  first  step  in 
the  potential  to  reduce  and  model  the  variability  in  materials  fati¬ 
gue  performance.  By  focusing  on  the  mechanics  and  physics  that 
control  the  short-life  tail  of  a  material’s  fatigue  capability,  it  ap¬ 
pears  possible  to  characterize  individual  contributions  to  life  limits 
with  a  high  level  of  fidelity,  potentially  leading  to  increases  in 
capabilities  to  control  and  fully  utilize  inherent  materials  proper¬ 
ties.  Two  primary  avenues  for  further  development,  maturation, 
and  transition  of  the  concepts  of  life-limiting  fatigue  are:  (i) 
through  the  integrated  design  and  certification  of  materials  and 
structures  for  future  systems  and  (ii)  through  improvements  to 
methods  of  life-cycle  management  and  prognosis. 

As  discussed  above,  beginning  with  the  residual-stress-free 
case,  one  can  start  to  isolate  effects  of  microstructure  on  inherent 
fatigue  variability.  This  special  case  presents  the  opportunity  to 
identify  life-limiting  microstructural  neighborhoods  and  effects 
of  microstructural  gradients  on  fatigue.  As  shown,  the  inherent  life 
limits  under  low  cycle  fatigue  appear  to  be  governed  by  the  growth 
of  a  dominant  crack  from  an  initial  crack  size  given  by  the  micro- 
structural  feature  that  controls  crack  formation.  Thus,  the  approach 
of  segmenting  total  lifetime  into  the  components  of  initiation, 
small-,  and  large-crack  growth  offers  a  very  systematic  process 
for  characterizing  and  modeling  the  many  additional  effects  that 
are  often  difficult  to  quantify  independently.  Such  models  should 
provide  new  interpretations  of  fatigue  behavior  that  have  tradi¬ 
tionally  been  based  primarily  on  extrapolations  from  mean-fatigue 
properties.  This  could  be  significant  in  some  instances,  since  min¬ 
imum  and  mean  fatigue  responses  differ  based  on  microstructural 
and  loading  variables. 

In  light  of  the  rapidly  growing  science  and  technology  offered 
by  Integrated  Computational  Materials  Engineering,  the  modeling 
of  life-limiting  fatigue  of  materials  presents  many  opportunities 
for  the  development,  transition,  and  implementation  of  new 
understanding  for  both  military  and  commercial  applications, 
including: 

•  Quantitative  2-  and  3-dimensional  characterization  and  repre¬ 
sentation  of  life-limiting  microstructural  features  vs.  nominal 
and  off-nominal  microstructures,  including  characterization  of 
microstructural  hierarchies  that  control  failure. 

•  Nondestructive  characterization  of  materials  deformation,  dam¬ 
age  states,  microstructural  features,  and  residual  stresses. 
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•  High  resolution  microscale  testing  of  microstructural  features, 
including  digital  image  correlation  for  precision  characteriza¬ 
tion  of  crystallographic  mechanisms  of  deformation  and 
damage. 

•  Crystal  plasticity  modeling  of  microstructural  features,  includ¬ 
ing  life-limiting  local  neighborhoods  and  gradients. 

•  Analysis  and  modeling  of  effects  of  surface  and  bulk  residual 
stresses. 

•  Effects  of  geometric  stress  concentrations,  such  as  bolt  holes 
and  notches. 

•  3-dimensional  fracture  mechanics  of  complex  geometric 
features. 

•  Effects  of  complex  mission  loading  and  histories,  including  ther¬ 
momechanical  fatigue  and  dwell  hold  times. 

•  Advanced  validation  materials  science,  including  selective  anal¬ 
ysis  of  the  damage  states  of  materials  from  retired  components 
to  provide  model  feedback  and  validation. 

•  Information  integration  and  data  mining. 

•  Methods  for  modeling  uncertainty  and  confidence  bounds  in 
minimum  fatigue  life  due  to  limited  availability  of  data  [67] 
and  for  statistical  optimization  of  testing. 

•  Methods  of  material  state  awareness  supporting  condition 
based  management  +  prognosis. 

•  Probabilistic  modeling  of  multi-scale  domains:  microstructure, 
geometric  feature,  component,  and  system. 

•  Model-based  probabilistic  simulation  and  validation  of  materi¬ 
als  and  components  designs. 

5.  Conclusions 

This  paper  used  electropolished  specimens  of  the  high-strength 
a  +  p  titanium  alloy  Ti-6Al-2Sn-4Zr-6Mo  to  explore  fundamental 
variabilities  and  uncertainties  in  microstructurally  based  fatigue 
life  limits  under  stresses  and  lifetimes  representative  of  engine  ro¬ 
tor  design  (nominally  8000  major  cycles,  (TACs)).  In  this  life  range, 
it  is  well  known  that  the  total  fatigue  lifetime  (NT)  may  be  parti¬ 
tioned  as  (Eq.  (2)): 

Nj  =  N\  Nsc  +  Nlc 

where  Nit  Nsc,  and  NLC  are,  respectively,  the  numbers  of  cycles  spent 
in  crack  initiation,  small-crack  growth,  and  large-crack  growth  to 
failure.  In  the  test  results  presented  here,  N[  was  often  the  dominant 
term  contributing  to  the  mean  fatigue  lifetime,  but  the  mean  fatigue 
behavior  sometimes  also  encompassed  a  shorter-life  fatigue  re¬ 
sponse  wherein  N[  was  small  (approached  zero).  This  minimum 
vs.  mean  bimodal  fatigue  behavior  appeared  to  respond  differently 
to  stress  level,  which  may  be  expected  depending  on  the  relative 
contribution  of  N[  to  the  total  fatigue  lifetime.  In  related  research, 
we  have  also  observed  that  the  trends  in  minimum  vs.  mean  fatigue 
response  may  differ  as  a  function  of  microstructural  and  usage 
variables. 

For  turbine-engine  rotor  materials,  traditional  design  practice 
defines  the  safe-life  fatigue  limit  as  a  statistical  minimum  life  to 
form  a  small  surface  or  corner  crack.  Alternatively,  from  a  materi- 
als-microstructure  perspective,  one  might  define  the  minimum 
size  of  a  small  crack  as  the  dimension  of  the  microstructural  fea¬ 
tures  in  which  cracks  initiate,  which  for  the  Ti-6246  alloy  are  pri¬ 
mary  alpha  particles  approximately  4  pm  in  diameter.  Under  the 
latter  definition  of  cracking,  our  consistent  results  on  repeated 
tests  indicated  that  the  minimum  of  N[  approached  0  cycles,  which 
meant  that  the  minimum  lifetime  was  spent  almost  completely  in 
the  growth  of  a  crack  that  begins  on  the  microstructural  scale. 
Although  the  probability  of  occurrence  of  N[  -►  0  may  be  reduced 
by  factors  such  as  the  imposition  of  surface  residual  stresses 
through  shot  peening,  in  the  stress  range  of  rotor  design,  it  appears 


difficult  to  preclude  the  rare  occurrence  of  secondary  damage, 
which  also  leads  to  N[  -►  0. 

Our  research  on  Ti-6246  and  other  high  performance  materials 
suggests  that  the  presumption  of  N[  -►  0  for  rotor  design  stresses 
may  be  of  general  applicability.  If  this  is  correct,  it  may  be  attrac¬ 
tive  to  reinterpret  conventional  approaches  and  strategies  for 
life-cycle  design  and  management  of  turbine-engine  rotor  materi¬ 
als  and  components. 
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